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G^ : Abstract 

I ■ The total hadronic cross-section a-y^{W) for the interaction of real photons, 77 — > hadrons, 

^. is measured for 77 centre-of-mass energies 10 < W^ < 110 GeV. The cross-section is extracted 

r-| ! from a measurement of the process e"'"e~ — > e"'"e~7*7* — > e"'"e~ + hadrons, using a luminosity 

^ I function for the photon flux together with form factors for extrapolating to real photons [Q"^ = 

Sp ' ^ GeV^). The data were taken with the OPAL detector at LEP at e"'"e~ centre-of-mass energies 

Jh ' y/^ee ~ 161, 172 and 183 GeV. The cross-section ajy{W) is compared with Regge factorisation 

. . . . and with the energy dependence observed in 7p and pp interactions. The data are also compared 

to models which predict a faster rise of a^^{W) compared to 7p and pp interactions due to 

additional hard 77 interactions not present in hadronic collisions. 
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1 Introduction 



At high 77 centre-of-mass energies W = y/s , the total hadronic cross-section a^^ for the 
production of hadrons in the interaction of two real photons is expected to be dominated by 
interactions where the photons have fluctuated into a hadronic state. Measuring the a/s 
dependence of a.y^ should therefore improve our understanding of the hadronic nature of the 
photon and the universal high-energy behaviour of total hadronic cross-sections. 

Before data from LEP became available, the total hadronic 77 cross-section had only 
been measured for 77 centre-of-mass energies W below 20 GeV by PLUTO |T]|, TPC/27 |0], 
PEP/27 ^^^d the MDl experiment 0], in a kinematic region where the expected high-energy 
rise of the total cross-section could not have been observed. Using LEP data taken at e"'"e~ 
centre-of-mass energies y^^^ = 130 — 161 GeV, L3 p has demonstrated that the total hadronic 
77 cross-section in the range 5 <W < 75 GeV is consistent with the universal Regge behaviour 
of total cross-sections. 

Processes with a pointlike coupling of the photon to quarks are absent in hadron-hadron 
collisions. This additional hard component in photon interactions is therefore expected to 
lead to a different energy dependence of the total cross-section for photon-induced interactions 
in comparison to hadron-hadron scattering. Models |^, |^ based on perturbative Quantum 
Chromodynamics (pQCD) encompass this by including additional photon interactions, usually 
denoted "direct" and "anomalous", in addition to the interactions which are described by 
the Vector Meson Dominance model (VMD). In VMD models, the photon fluctuates into a 
bound state vector meson. In Regge models, a different energy dependence of the cross-section 
can be obtained either by universality breaking effects or by introducing an additional hard 
pomeron [§ . Thus far no sign of such a different energy dependence has been experimentally 
established in comparison of the total 7p cross-section measured by the HERA experiments |^, 
10|] with total pp cross-sections. Any such effect would be expected to be more pronounced in 



77 interactions, since here one has two photons in the initial state. 

In this paper, we present a measurement of the total hadronic 77 cross-section in the 
range 10 < VT < 110 GeV using data taken by the OPAL detector at LEP at -/Sgg = 161, 



172 and 183 GeV. The integrated luminosities are 9.9, 10.0 and 54.4 pb~^, respectively. At 
these energies above the Z° resonance, hadron production is dominated by photon-photon 
collisions and background from other processes, e.g. e+e~ annihilation, is expected to be small. 
The photon-photon events are selected by a series of cuts intended to exclude backgrounds, 
especially from the e"'"e~ annihilation and 77 —>■ i~^i~ channels [i G {e,/i, r}). In addition, 
an anti-tagging condition is applied, requiring that no scattered electron^ was detected. Most 
of the photons therefore carry only a small negative four-momentum squared, Q^, and can be 
considered to be quasi-real {Q^ ~ GeV^). The differential cross-section da/dW is measured 
for the process e"'"e~ — >■ e"'"e~ + hadrons, where W is the invariant mass of the hadronic 
system. From this cross-section the total hadronic cross-section a^^{W) for the interaction of 
real photons, 77 -^ hadrons, is extracted as a function of W, using a luminosity function for 
the photon flux and form factors for the extrapolation to Q"^ = GeV^. 



2 The OPAL detector 



A detailed description of the OPAL detector can be found in Ref. [11], and therefore only a 



brief account of the main features relevant to the present analysis will be given here. 

The central tracking system, covering the polar angle range | cos6\ < 0.73, is located inside 
a solenoidal magnet which provides a uniform magnetic field of 0.435 T along the beam axis|. 
The magnet is surrounded in the barrel region (| cos6'| < 0.82) by a lead glass electromagnetic 
calorimeter (ECAL) and a hadronic sampling calorimeter (HCAL). Outside the HCAL, the 
detector is surrounded by muon chambers. There are similar layers of detectors in the endcaps 
(0.82 < I cos^l < 0.98). The small-angle region from 47 to 140 mrad around the beam pipe on 
both sides of the interaction point is covered by the forward detectors (FD) and the region from 
25 to 59 mrad by the silicon-tungsten luminometers (SW). From 1996 onwards, relevant to the 
data presented in this paper, the lower boundary of the SW acceptance has been increased to 
33 mrad following the installation of a low-angle shield to protect the central detector against 
synchrotron radiation due to the increased LEP e"'"e~ beam energies. 

Starting with the innermost components, the tracking system consists of a high precision 
silicon microvertex detector, a vertex drift chamber, a large volume jet chamber with 159 
layers of axial anode wires and a set of z chambers measuring the track coordinates along 
the beam direction. The transverse momenta px of tracks are measured with a precision 



parametrised by (Jp^/p'-£ = JO.02'^ + (0.0015 ■ Pt^ {pt in GeV/c) in the central region. In this 
paper, "transverse" is always defined with respect to the z axis. The jet chamber also provides 
measurements of track energy loss, dE/dx, which are used for particle identification |[Tl[ . 



The barrel and endcap sections of the ECAL are both constructed from lead-glass blocks, 
with a depth of 24.6 radiation lengths in the barrel region and more than 22 radiation lengths 
in the endcaps. The FD consist of cylindrical lead-scintillator calorimeters with a depth of 24 
radiation lengths divided azimuthally into 16 segments. The electromagnetic energy resolution 



^Positrons are also referred to as electrons 

^In the OPAL coordinate system the z axis points in the direction of the e~ beam. The polar angle 9, the 



azimuthal angle (f> and the radius r denote the usual spherical coordinates. 



is about 18%/\/E, where E is the energy in GeV. The SW detectors [12] each consist of 19 
layers of silicon interleaved with 18 layers of tungsten, corresponding to a total of 22 radiation 
lengths. Each silicon layer consists of 16 wedge shaped silicon detectors. The electromagnetic 
energy resolution is about 25% /^/E {E in GeV). 



3 Kinematics 



A schematic diagram of the two-photon process is shown in Fig. |l|. The kinematics of the 
process e"'"e~ — >■ e"'"e~ + hadrons at a given -y/igg can be described by the negative square of the 
four-momentum transfers, Qf = —qf, carried by the two {i = 1, 2) incoming virtual photons 
(7*) and by the square of the invariant mass of the hadronic final state, W"^ = s^^ = (gi + 52)^- 
The four-momenta of the electrons before and after the interaction are denoted by Pi and p'^, 
respectively. Each Qf is related to the electron scattering angle 6'^ relative to the beam direction 
by 

Ql = -fe - P-? ^ 2E,E:(1 - cos^O, (1) 

where Ei and E'^ are the energies of the beam electron and the scattered electron, respectively. 
Events are only included in the analysis if they do not contain scattered electrons (either single- 
tagged or double-tagged events). This anti-tagging condition defines an effective upper limit 
on the values of Qf for both photons. This condition is either met if the scattering angle 6'^ of 
the electron is less than 33 mrad, defined by the angle between the beam axis and the inner 
edge of the acceptance of the SW detector, or if the energy of the scattered electron is smaller 
than the minimum energy of 20 GeV required for the tagged electron in SW or 40 GeV in FD. 



4 Event selection 

Two-photon events are selected with the following set of cuts: 

• The visible invariant mass calculated from the position and the energy of the clusters 
measured in the ECAL has to be greater than 3 GeV. 

• The sum of all energy deposits in the ECAL and the HCAL has to be less than 45 GeV 
in order to reject e"'"e~ annihilation events. 

• At least 2 tracks must have been found in the tracking chambers. A track is required 
to have a minimum transverse momentum of 120 MeV/c, at least 20 hits in the central 
jet chamber, and the innermost hit of the track must be inside a radius of 60 cm with 
respect to the z axis. The point of closest approach to the origin in the r0 plane must 
be less than 20 cm in the z direction and less than 1 cm in the rep plane. Tracks with 
a momentum error larger than the momentum itself are rejected if they have less than 
80 hits. The number of measured hits in the jet chamber must be more than half of the 
number of possible hits, where the number of possible hits is calculated from the polar 
angle 6 of the track, assuming that the track has no curvature. 



• The transverse momentum of the event measured in the ECAL and the FD has to be less 
than 5 GeV/c. 

• No track in the event has a momentum greater than 30 GeV/c. 

• To remove events with scattered electrons in the FD or in the SW calorimeters, the total 
energy sum measured in the FD has to be less than 40 GeV and the total energy sum 
measured in the SW less than 20 GeV. This cut also reduces the contamination from 
multihadronic e^e~ annihilation events with their thrust axis close to the beam direction. 

• The background due to beam-gas or beam-wall interactions is reduced by the following 
requirements. The radial distance of the primary vertex from the beam axis has to be 
less than 3 cm. To estimate the z position of the primary vertex for photon-photon 
events with typically low multiplicity, we calculate the error-weighted average (zq) of 
the z coordinates of all tracks at the point of closest approach to the origin in the r0 
plane. The background due to beam-gas or beam-wall interactions is further reduced by 
requiring |(-2o)| < 10 cm and that the net charge of an event, calculated by adding the 
charges of all tracks, is less or equal three. 

For the remaining events we determine the visible energy i?vis and the longitudinal com- 
ponent, Pl, and the transverse component, Pt? of the momentum vector of the hadronic final 
state. These quantities are calculated after a matching algorithm is applied to the data, in order 
to avoid double-counting of particle momenta. The matching algorithm uses all the information 
from the ECAL and the HCAL, the FD and the SW calorimeters, as well as from the tracking 
system. If a calorimeter energy cluster is associated to a track, the cluster energy is compared 
to the expected energy response /(p) of the calorimeters for the track with momentum p. To 
calculate the energy associated to a track, the pion mass is assumed. The cluster is rejected if 
the energy of the cluster is less than expected from the track energy. If the cluster energy E 
exceeds the expected energy by more than what is expected from the resolution, the energy of 
the cluster is reduced to E — f{p). In this case the track momentum and the reduced energy of 
the cluster are taken separately. The output of the matching algorithm is an array of energies 
and momenta (i?h,Ph) which are used to calculate the visible invariant mass Wyis'. 

^™ = [T.E^J-[T.P^^' (2) 

= El,-Pl-Pl (3) 

A cut |PL/i?vis| < 0.85 is applied, since beam-gas or beam-wall events tend to accumulate at 
high values of IPlZ-EvIsI- 

In order to reject events with only leptons in the final state, additional requirements have 
to be fulfilled for events with rict = 2, where rich is the number of tracks. Each of the two 
tracks must have at least 20 dE/dx hits in the central jet chamber and the dE/dx probability 
must be smaller than 10% for the electron and for the muon hypothesis. The thrust of the 
event, calculated in the laboratory system from the output of the matching algorithm, has to 
be below 0.98. 

We use data corresponding to an integrated luminosity of 74.4 pb^^ collected with the OPAL 
detector in the years 1996 and 1997. The integrated luminosity is 9.9 pb~^ at v/iee ~ l^l GeV, 
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10.0 pb ^ at i/igg = 172 GeV and 54.4 pb ^ at i/See ~ ^^^ GeV. The error on the luminosity 
is less than 1%. After applying the above cuts, 23250 events remain at \/s^^ = 161 GeV, 25643 
at ^^^ = 172 GeV and 144147 at ^/s^^ = 183 GeV. 



5 Monte Carlo simulation 



The Monte Carlo generators PYTHIA 5.722 and PHOJET 1.05c are used to simu- 
late photon-photon interactions with Qf and Q2 < 4.5 GeV^. The photon-photon generator 
PYTHIA is based on a model by Schuler and Sjostrand [0 and PHOJET has been developed 
by Engel based on the Dual Parton model (DPM) [|I4|. Both generators simulate the pro- 
cess (e"'"e~ — >• e"'"e~ + hadrons) in two stages, firstly e"'"e~ — >• e"'"e~77 and then 77 -^ hadrons. 
The probability of the beam electron emitting a photon is modelled by the Equivalent Photon 
Approximation (EPA) [0. 



Soft processes like quasi-elastic scattering (77 -^ VV, where l^ is a vector meson), single- 
diffractive scattering (77 -^ VX, where X is a low mass hadronic system) or double-diffractive 
scattering (77 — >■ X1X2) are modelled by both generators. The cross-sections are obtained by 
fitting a Regge parametrisation to pp, pp and 7p data and by assuming Regge factorisation, 
i.e. universal couplings of the pomeron to the hadronic fluctuations of the photon. In both 
generators the quasi-elastic cross-section is about 5 — 6%, the single- diffract ive cross-section 
about 8 — 12% and the double-diffractive cross-section about 3 — 4% of a^^ for W > 10 GeV. 

The transition from soft to hard interactions is defined by the transverse momentum of 
the primary produced partons. For the hard interactions it is assumed that the cross-section 
can be factorized into parton distribution functions which give the probability to find a parton 
(quark, gluon) in the photon and matrix elements for the hard subprocess. All possible hard 
interactions of quarks, gluons and photons are simulated using leading order (LO) matrix 
elements. As default the SaS-lD parametrisation of the parton distribution functions |16[ is 
used in PYTHIA and the LO GRV parametrisation |T7[ in PHOJET. 

The fragmentation and decay of the parton final state is handled in both generators by the 
routines of JETSET 7.408 P]. Initial- and final-state parton radiation is included in the leading 
logarithm approximation. Both generators include multiple interactions of the remnants of the 
initial photons. 

The two-photon mode of PYTHIA simulates the interactions of real photons with Q\.,Q\ = 
GeV^. The virtuality of the photons, defined by Q^, enters only through the EPA in the 
generation of the photon energy spectrum, but the electrons are scattered at zero angle. In 
PHOJET the Q^ suppression of the total 77 cross-section is parametrised using Generalised 
Vector Meson Dominance (GVMD). The Q^ dependence of the quark and gluon densities of the 
virtual photon and additional Q^ dependent suppression factors for diffractiveQ processes are 
also taken into account [^ . The Q^ dependent transverse momenta of the scattered electrons 
are also simulated. 

All signal and background Monte Carlo samples are generated with full simulation of the 



■^In this paper we refer to the sum of quasi-elastic, single- and double-diffractive events as diffractive events 
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OPAL detector [18|. They are analysed using the same reconstruction algorithms as are applied 
to the data. The background from e+e" annihilation events e+e" -^ (7/Z°)* -^ qq(7) is gen- 
erated with PYTHIA M. The leptonic two-photon background processes e+e~ -^ e+e'r+r", 



e"'"e~ — >■ e"'"e^/x+/i^ and e^e~ —>■ e"'"e~e+e~ are simulated with VERMASEREN |]19[. The contri- 
bution from other background processes is negligible. The Monte Carlo simulated background 
is less than 1.6% of the total number of selected events . This does not include beam-gas and 
beam-wall interactions which are estimated to contribute about 2% to the total number of 
selected events. 



Deep-inelastic e^y (= 7*7) events are generated with HERWIG 5.9 pOl. From the Monte 
Carlo it is estimated that after all cuts about 1.5% of the remaining events are e7 processes 
with Taa:x{Qf,Q2} > 4.5 GeV^ for ^/s^^ = 183 GeV. The rate of events where both Qf are 
larger than 4.5 GeV^ is negligible. 



6 Unfolding of the hadronic cross-section 

In a first step, the differential cross-section daee/dW for the process (e"''e~ -^ e'^e~ + hadrons) 
is obtained from the W^is distribution. The measured V^vis distribution is shown in Fig. |^ for 
the data taken at i/igg = 183 GeV. The distribution, which falls smoothly over 5 orders of 
magnitude, is well described by the Monte Carlo simulations which have been normalized to 
the number of data events after subtracting the Monte Carlo expectation for background and 
e7 events. 

The selection efficiency, defined as the ratio of the number A^gei of selected Monte Carlo 
events to the number iVgen of generated events at a given generated W, is shown in Fig. ^ It 
rises from about 10% at PF = 10 GeV to an almost constant plateau of about 65% for the 
PHOJET events with W > 40 GeV and it decreases again for very high W. The selection 
efficiency for the PYTHIA events is about 15% lower at W^ = 40 GeV and it approaches the 
PHOJET selection efficiency at high W. The selection efficiency for diffractive events simulated 
with PHOJET is much higher than for the diffractive events simulated with PYTHIA. 

The relation between W^is and the generated W for all selected PHOJET and PYTHIA 
Monte Carlo events is shown in Fig. ^. The finite resolution of the W measurement is given by 
the standard deviations of the l^vis distributions in each bin of W which are plotted as vertical 
bars. The W resolution is consistent for both generators, but differences for the average VFvis 
as a function of the generated W are observed at high W. The main energy losses are caused 
by hadrons which are emitted at small polar angles 6; they are either lost in the beam pipe, or 
they are only detected with low efficiency in the electromagnetic calorimeters in the forward 
regions (FD and SW). 

The background determined by the Monte Carlo is first subtracted from the data. Then, 
the unfolding of the resolution effects, as well as the correction for the detector acceptance 
and the selection cuts are done with the program GURU |21||. This program for regularised 
unfolding is based on the Singular Value Decomposition (SVD) method. For systematic checks. 



the unfolding program RUN |22] has also been used. 



The differential cross-sections for the three beam energies, daec/dW, after unfolding are 
given in Fig. |^ and Table Q. Bin-to-bin correlations are sizeable, since the chosen bin size is not 
much larger than the resolution. The size of the correlations also depends on the regular isat ion 
procedure of the unfolding. The covariance matrix obtained from the unfolding is given in 
Table |. 

The differential cross-section dage/dW of the process (e+e^ -^ e+e^ + hadrons) can be 
translated into the cross-section a^^ for the process (77 — >■ hadrons) using the luminosity func- 
tion L^^ for the photon flux [|l^. The cross-section for real photons is derived by using form 
factors -F(Q^) which describe the Q^ dependence of the hadronic cross-section. In every W bin 
AWi we determine 

(4) 
where yi and 1/2 denote the fraction of the beam energy carried by the photons with yiy2 ~ 
VT^/see (neglecting Q\ and Q"^)- The cross-section a^^ is given at the bin centre, since the 
deviation of Wl from the bin centre due to the finite bin width is found to be small. 

The luminosity function L^^ and the form factors F{Q'^) for the various W bins are obtained 
from the program PHOLUM [^ which performs a numerical integration for each W bin over the 
unmeasured phase space (Qi, Q\ and y ranges). PHOLUM takes into account both transverse 
and longitudinally polarized photons. The form factors are used in the GVMD approximation 
of Ref. p3[. The difference between the extrapolation to Qi = and Q2 = is about 7% of 
a^^ if the GVMD model is compared to a simple p° form factor |]^. This uncertainty is not 
included in the systematic error of the measurement. 

In the analysis, the e7 events simulated by HERWIG are subtracted from the data and 
the photon flux is therefore calculated with a cut on maxjQf, Q^}- I^ order to check this 
procedure, the analysis was also performed without subtracting the Monte Carlo e7 events. In 
this case, the photon flux has to be calculated without a cut on max{Q^, Q\}. The uncertainty 
is estimated by comparing these procedures and by using PYTHIA instead of HERWIG for 
the modelling of the e7 background. The resulting uncertainty on a^^ is about 1% and it is 
therefore neglected. 

Assuming the Q\, Q\ and W dependence of the total hadronic cross-section for virtual 
photons, 0"^.^. (Q^, Q\^ VT), to factorize, based on a simple GVDM ansatz, the W dependence is 
preserved when extrapolationg to a^^{yV\ Although there are some events in the tails towards 
higher Q\ or Q\ in the 77 data of this experiment, extending to several GeV^, the bulk of the 
data are at very low Q\ and Q\^ so a significant contribution from the tail is unlikely. The 
medians of the Q\ and Q\ distributions are of the order 10~^ GeV^ (taken from Monte Carlo). 

Radiative corrections like multiple photon emission off the incoming electrons are not in- 



cluded in the Monte Carlo generators. They are expected to be small |2J] and the effect of the 
radiative corrections should be much reduced by using the hadronic final state to calculate the 
kinematics, i.e. the hadronic invariant mass W ^ of the event. 

The three data samples at a/s^p = 161, 172 and 183 GeV were independently analysed and 
the results for the total hadronic two-photon cross-section o^^ are found to be in agreement 
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within 1-2 standard deviations of the statistical error. Furthermore, no systematic trend in the 
W dependence of cr^^iW) is observed as a function \/Sgg. The total cross-sections are therefore 
averaged using as weight the corresponding integrated luminosities (Tab. |^). 



7 Systematic errors 

Several distributions of the data are compared to PYTHIA and PHOJET after detector sim- 
ulation in order to study whether the general description of the data by the Monte Carlo is 
sufficient to use the Monte Carlo for the unfolding of the cross-section. The Monte Carlo distri- 
butions are all normalized to the number of data events after the Monte Carlo expectation for 
background and ej events were subtracted from the data. Without this normalisation, using 
the cross-section predicted by the Monte Carlo generators, the number of selected events is 
about 10% smaller than in the data for PHOJET and about 10% larger than in the data for 
PYTHIA. 

In both Monte Carlo models about 20% of the cross-section is due to diffractive events in 
which the final state hadrons go strongly forward or backward into those parts of the detector 
which have the smallest acceptance. This fraction is almost independent of W for W > 10 GeV. 
The selection efficiency for the diffractive events is small and, although the generated rate is 
almost the same in both models, different modelling of the diffractive events leads to very 
different selection efficiencies. For a W = 70 GeV only about 6% of all generated diffractive 
events are selected in PYTHIA, whereas about 20% are selected in PHOJET (Fig. 0). The 
detector correction therefore has to rely heavily on the Monte Carlo simulation for this class of 
events. 

In order to study the modelling of the diffractive events using the data, we have plotted 
the maximum rapidity gap Ar^max between the pseudorapidities rj = —In tan 9/2 of any two 
particles, neutral or charged, found by the matching algorithm in Fig. |^. Diffractive events are 
expected to have larger Ar^max due to the colour- singlet exchange |^. The data are compared 
to the PHOJET and PYTHIA simulations. Both models underestimate the Ar^max distribution 
at large Ay^max with PHOJET being closer to the data than PYTHIA. It was checked that the 
transverse momenta of the scattered electrons, which are simulated in PHOJET but not in 
PYTHIA, have only a small effect on this distribution. 

Significant discrepancies are also found in the distribution of the charged multiplicity n^h 
(Fig. 1^) and the distribution of the thrust variable, T, (Fig. H). Both Monte Carlo models 
significantly underestimate the fraction of low-multiplicity events [rich < 6) and overestimate 
the fraction of high-multiplicity events in comparison to the data and there are also more events 
with large thrust (T > 0.925), in the data. It should be noted that increasing the fraction of 
diffractive events by factors of two or more in the Monte Carlo does not lead to a significant 
improvement in these comparisons. 

The energy -Esw measured in the silicon-tungsten luminometers (SW) is shown in Fig. |^ for 
all selected events with Esw > 1 GeV and the energy EpD measured in the forward detectors 
(FD) is shown in Fig. |1^ for all selected events with E'fd > 2 GeV. At low i^sw both Monte 
Carlo models lie above the data, but the reasonable agreement of data and Monte Carlo at 
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large -Esw and Ep^ shows that the remaining background from multihadronic e^e~ annihilation 
events and deep-inelastic e7 events is small and that this remaining background is reasonably 
well described by the Monte Carlo. This implies that there are also no events left with an 
off-momentum beam electron which was scattered upstream hitting SW or FD. 

Finally, we plot the ratios PtZ-^^vIs and PlZ-E'vIs of the transverse and longitudinal com- 
ponents of the momentum vector of the hadronic system to the visible total energy E, 



VIS 



(Figs. |TI|,|12D. Data and Monte Carlo are in good agreement. The small number of events 
at large PxZ-E'vis are expected to be mainly due to background processes. Studies of beam-gas 
and beam-wall events show that most of these events have |PL/-Evis| > 0.85. The background 
conditions were different at the three beam energies. In the data taken at ^/s^^ = 183 GeV the 
beam related background peaked mainly at positive PlZ-EvIs > 0.85. This explains the small 



asymmetry in Fig. |T2]c 



Based on these observations, the following systematic errors are taken into account in the 
measurement of the cross-sections for every beam energy, separately (Tab. |lD: 

• In most of the distributions, both Monte Carlo models describe the data equally well and 
there is no reason for preferring one model over the other for the unfolding of the data. We 
therefore average the results of the unfolding. The difference between this cross-section 
and the results obtained by using PYTHIA or PHOJET alone are taken as the systematic 
error due to the Monte Carlo model dependence of the unfolding. 

• An additional error due to the uncertainties of the modelling of the diffractive processes in 
the Monte Carlo is taken into account. Since there is large uncertainty on the diffractive 
77 cross-section derived from the HERA measurements |^, |I^, |2B|, we have increased the 



percentage of diffractive events from 18% to 27% in PHOJET which leads to an increase 
of a^j by 6%. Increasing the selection efficiency for diffractive events by a factor 2 leads 
to a decrease of a^-y by 6%. These variations of ±6% are used as systematic error. 

The uncertainty in the ECAL energy scale was estimated to be ±3% by comparing the 
energy distribution reconstructed in the ECAL for e"'"e~ annihilation events at i/igg = 
183 GeV with the Monte Carlo simulation. ECAL clusters of more than 10 GeV were 
excluded from this comparison in order to have a distribution of the energy per cluster 
which is similar to 77 events. The systematic error on the total cross-section was then 
estimated by varying the reconstructed ECAL energy in the Monte Carlo by ±3%. 

The electromagnetic calorimeters in the forward direction, SW and FD, are used in the 
Wvis measurement. A possible uncertainty in the energy scale and the detector simulation 
for hadrons reconstructed in SW or FD was studied by calculating and unfolding W^is 
without SW and FD information, respectively. The difference between a.y^{W) obtained 
without SW or FD information and a^^iW) obtained with the full detector is taken as 
the systematic error. 

The trigger efficiency was studied using data samples which were obtained using nearly 
independent sets of triggers. The trigger efficiency is defined as the ratio of the number 
of triggered and selected events to the number of selected events. On average, the trigger 
efficiency for the low W range, 10 < W < 35 GeV, is greater than 96% and it approaches 
100% for larger values of W. Only lower limits on the trigger efficiency can be determined 
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with this method and therefore no correction factor is apphed. However, the lower hmit 
on the trigger efficiency is taken into account as an additional systematic error. 

Studying vertex and net charge distributions, it is estimated that about 2% of the selected 
events could be due to beam-gas or beam-wall interactions. Hadronic photon-photon 
events, however, in coincidence with an off-momentum beam electron which was scattered 
upstream and hitting SW or FD are rejected by the SW and FD energy cuts. The fraction 
of photon-photon events rejected due to these coincidences is estimated to be less than 
2%. Taking into account both effects, a value of 2% is therefore taken as additional 
systematic error. 



The program GURU |£l| has been used for unfolding the W^is distribution. Since the 
distribution of the charged multiplicity rich is not well described by the Monte Carlo 
models, we have studied the influence of this discrepancy by unfolding the two-dimensional 
{Wvis, rich) distribution. No significant difference from the one- dimensional unfolding using 
GURU is found. The unfolding program RUN p2| can only be used for one- dimensional 



unfolding. The two unfolding methods implemented in RUN and GURU applied to the 
M^vis distribution also yield consistent results. Therefore no additional systematic error 
has been taken into account. 

• Since the background rate taken from Monte Carlo is only about 1.6%, a possible sys- 
tematic error is neglected. 

• The overall normalisation error due to the uncertainty on the luminosity measurement is 
less than 1% and is therefore also neglected. 

The different systematic errors are summed up in quadrature to obtain the total systematic 
error. For the total error, the statistical and the total systematic error are added in quadrature. 
The luminosity- weighted average values of the total cross-section a^-y{W) and errors for the 
different W bins are given in Table 0. 



8 Results and model comparisons 



The total cross-section for the process 77 -^ hadrons, a^^{W), is shown in Fig. 13 in the range 
10 < VT < 1 10 GeV. In the region W < 20 GeV, the OPAL measurement is consistent with the 
results from PLUTO [0], TPC/27 and PEP/27 within the large spread and experimental 
errors of these measurements. 

The OPAL measurements exhibit the rise in the W range 10 < VT < 110 GeV which is 
characteristic for hadronic cross-sections in this energy range. A similar rise was first observed 
by the L3 experiment 0, but their values of a^^ are about 20% lower than the OPAL measure- 
ment. L3 used PHOJET only for the unfolding, whereas for the OPAL measurement presented 
here the unfolding results of PHOJET and PYTHIA are averaged. The OPAL result obtained 
using only PHOJET is about 5 — 10% lower than the averaged result. 

Several models have been proposed to describe the energy dependence of hadronic cross- 
sections. One of the interesting questions for hadronic interactions of real photons is whether 
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they behave the same as hadrons or whether the additional hard contributions to the total cross- 
sections of photon-induced interactions lead to a faster rise of the total 77 and 7p cross-sections 
as a function of energy. Hence we performed a detailed study of the data in the framework of 
present models. 

We study the data within the framework of Regge theory. The total cross-sections for 



hadron-hadron and photon-proton collisions have been found to be well described [^, ^ by a 
Regge parametrisation of the form 

(Tab = XiABs'^ + YiABs"'' - Y^abs-"^', (5) 

where A and B denote the interacting particles and the centre-of-mass energy squared, s, is 
taken in units of GeV^. The first term in the equation is due to soft pomeron exchange and the 



other terms are due to C-even and C-odd reggeon exchange, respectively [|28[. The exponents 
ei, rji and 172 are assumed to be universal, whereas the coefficients Xiab and Fjab are process 
dependent. The values of the exponents were determined in Ref. [^ by a fit to the pp, pp, 
TT^p, K^p, 7p and 77 total cross-sections: 

ei = 0.095 ±0.002, r/i = 0.34 ± 0.02, and 772 = 0.55 ± 0.02. (6) 

The non-zero value of the exponent ei = 0.095 ± 0.002 predicts a slow rise of the total cross- 
section with energy. The fit in Ref. [^ is dominated by the hadron-hadron data. For 7p 



and 77 collisions I2 = 0, i.e. Eq. ^ reduces to the original form proposed by Donnachie and 
Landshoff ||2^. In this combined fit the available 77 data, i.e. not including the OPAL data 
presented here, were fitted in the range W > 4 GeV, yielding Xi^^ = (156 ± 18) nb and 



Fi^^ = (320 ± 130) nb |28 



Assuming factorisation of the pomeron term Xiab, the total 77 cross-section can be related 
to the 7p and pp total cross-sections at centre-of-mass energies y/s = \/s = V^pp larger 
than about 10 GeV, where the pomeron trajectory should dominate: 

a,, ^ ^. (7) 

Most models for the high-energy behaviour of a^^ are based on this factorisation assumption 
for the soft part of the cross-section. In order to predict cr^^ via Eq. ^ the fit values of Xipp, 
Yipp, Y2pp, Xi^p and Fi^p are taken from Ref. ||2^ together with the exponents given in Eq. |^. 
This simple factorisation ansatz gives a reasonable description of a^^, but a faster increase of 
the cross-section than predicted by ei = 0.095 cannot be excluded, as shown in Fig. 



The errors on the data points are dominated by systematic errors which are highly corre- 
lated. Therefore, the subsequent fits were made using the statistical error, i.e. the covariance 
matrix given in Tab. ^ and the x^ values are only calculated with the statistical errors. If 
the fits are repeated using the correlation matrix from the unfolding and the total errors, the 
X^ values are reduced by about a factor 200 but the fit results are essentially unchanged. In 
order to take the correlations due to the systematical errors fully into account, the fits were 
repeated for each systematic error source, by shifting the values of the cross-sections at each 
energy by the corresponding amount of the systematic error. The final systematic error on 
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the fit parameters is then calculated as the square root of the quadratic sum of the differences 
between the shifted and unshifted fits to the data. 

In all subsequent fits to our data, we fix the reggeon term by using the values given in 
Ref. |2^, rji = 0.34, Yi^^ = 320 nb and F277 = 0, since we have no data at low W to constrain 
the fit in this region. We first check the universality of the exponent ei by fitting Eq. ^ to the 
data, leaving the exponent ei and the coupling Xi^^ free. The results of this fit (denoted by fit 
1) are given in Table |^ and shown in Fig. ^^ together with the OPAL data points. Only the 
diagonal elements of the covariance matrix from Tab. ^ are shown as error bars. From fit 1 we 
obtain 

ei = 0.101 ± 0.004(stat)+!;;J^9(sys) (8) 

which is in agreement with the value ei = 0.095 it 0.002 which describes the hadron-hadron and 
7P data. 

In Ref. 1^, a scheme is proposed for analysing hadron, real-photon and virtual-photon 
interactions. This scheme is still based on Regge phenomenology, but introduces an extra term 
which can be identified with an additional hard pomeron, which has an intercept significantly 
larger than one. This hard pomeron is assumed to be responsible for the fast rise of the virtual 
photon-proton 7*p cross-section. In this model, the cross-section is given by 

^AB = X,abs'' + X2ABS'' + Y,abs-''\ (9) 

In a second fit to our data (denoted by fit 2), we fixed ei and rji to the values in Eq. ^ and €2 
to the value 0.418 proposed in Ref. to study the significance of the term X2 in Eq. ^ We 
obtain 

X2^^ = (0.5 ± 0.2(stat)+}:^(sys)) nb. (10) 

Hence, we find that within the precision of our data this additional term is not required. 

Total cross-sections are also described by QCD inspired models. The total cross-sections 
obtained with some of the models are shown in Fig. |13|. All these models predict a steeper rise 
of a^^(W) than the simple factorisation ansatz based on the pp and 7p cross-sections which 
is not observed within the uncertainty of our data. Schuler and Sjostrand |13| give a total 
cross-section for the sum of all possible event classes in their model of 77 scattering where the 
photon has a direct, an anomalous and a Vector Meson Dominance (VMD) component. The 
direct and anomalous components lead to additional hard interactions which are calculated to 
leading order in pQCD. Schuler and Sjostrand consider the spread between this prediction and 
the simple factorisation ansatz as a conservative estimate of the theoretical band of uncertainty. 

We also plot the prediction of Engel and Ranft 0] which is implemented in PHO JET and an 
eikonalised mini-jet model by Godbole and Panchieri |^ which uses the GRV parton densities 
of the photon and a transverse momentum cut-off of 2 GeV/c for mini-jet production. The 
soft part of the cross-section is derived from 7p data. Another eikonalized mini-jet model, 
which assumes simple relations between photon and hadron-induced partonic cross-sections, 
also predicts cross- sect ions which are lower by about 20% than the data [pT| . 
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9 Conclusion 



We presented a measurement of the total hadronic cross-section a.yy{W) for the interaction of 
two real photons, 77 -^ hadrons, in the range 10 < W^ < 110 GeV for data taken at \^^^ = 161, 
172 and 183 GeV. The cross-section is in good agreement with a simple factorisation ansatz 
based on 7p and pp data. 

The energy dependence of the total cross-section has been studied using Regge model 
parametrisations. We observe the high-energy rise of the total cross-section which is typi- 
cal for hadronic interactions. We find that previous global fits to the y/s dependence of o"pp, 
(T-yp and a^^ give a good description of our results. The fit of the exponent in the soft pomeron 
term yields 0.101 ± 0.004+[]:[j?^ compared with 0.095 ± 0.002 from the global fit. Within the 
uncertainty of our measurement, no indication for a faster rise of the 77 cross-section than in 
hadron-hadron or 7p interactions is observed. 

Further improvements of the description of the hadronic final state by Monte Carlo models 
are necessary to reduce the systematic error of the measurement. It will also be important to 
gain a better understanding of the diffractive processes in 77 scattering. 
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W^-range [GeV] 


dcTec/diy [pb/GeV] 


da^JdW [pb/GeV] 


da^JdW [pb/GeV] 




v/i^^ = 161 GeV 


v/i^^ = 172 GeV 


v/i^^ = 183 GeV 


10- 20 


294.0 ± b.Gtlli 


319.4 ±5.5li| 


320.6 ± 2.4111 


20- 35 


89.0 ±1.91 II 


96.1 ± 1.9t\^oi 


104.4 ±o.9l|l:i 


35- 55 


33.5 ± 0.81 1? 


36.4 ±0.81 li 


39.4 ±0.41 f^g 


55- 80 


11.2 ±0.31 li 


13.9 ±0.41 li 


14.8 ±0.21 li 


80 -110 


3.5 ± 0.21 °:5^ 


4.8±0.2l ^:? 


5.3±0.ll[Ji 



Table 1: The differential cross-section do"ee/dH^ at v^,,^ = 161, 172 and 183 GeV for the 
anti-tagged two-photon events. The first error is statistical and the second systematic. 



W^-range [GeV] 


10- 20 


20- 35 


35- 55 


55- 80 


80-110 


10- 20 


6.75 


1.75 


-2.86 


-2.46 


0.84 


20- 35 


1.79 


8.72 


3.97 


-3.77 


-5.68 


35- 55 


-2.86 


3.97 


13.72 


10.07 


-6.19 


55- 80 


-2.46 


-3.77 


10.07 


26.85 


20.42 


80 - 110 


0.84 


-5.68 


-6.19 


20.42 


54.88 



Table 2: The covariance matrix of the statistical errors on cr^^{W) obtained from the unfolding 
The units are nb^. 
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W-r&nge [GeV] 


10- 20 


20- 35 


35- 55 


55- 80 


80- 110 


cr^-y [nb] 


362 


372 


414 


439 


464 


stat. error 


± 3 


± 3 


± 4 


± 5 


± 7 


MC model 


±21 


±30 


±29 


±29 


±53 


diffraction 


±22 


±22 


±25 


±26 


±28 


ECAL 


±16 


± 9 


± 9 


±10 


±14 


no FD 


- 1 


- 9 


- 5 


±13 


±41 


noSW 


- 4 


± 2 


± 7 


±10 


±11 


trigger 


±14 


± 9 


± 8 


± 7 


± 5 


beam-gas 


± 7 


± 7 


± 8 


± 9 


± 9 


total syst. 


+37 
-35 


+40 
-40 


+42 
-41 


+45 
-41 


+75 
-62 


total error 


+38 
-35 


+40 
-40 


+42 
-41 


+45 
-41 


+76 
-62 



Table 3: The total hadronic two-photon cross-section a^^ and the contributions from the various 
systematic errors (in nb) 





Xiyy [nb] 


ei 


X2^'y [nb] 


£2 


Yi.yy [nb] 


V 


xVndf 


fit 1 


180 ± 5li 


0.101 ± o.oo4t[!:°25 







320 


0.34 


68/3 


fit 2 


182 ± 3li 


0.095 


0.5±0.2li:^ 


0.418 


320 


0.34 


65/3 


PDG [H 


156 ±18 


0.095 ±0.002 







320 ± 130 


0.34 ±0.02 





Table 4: Results of the various fits of Regge type parametrisations to the total 77 cross-section. 
If no error is given, the parameter was fixed in the fit. The values of x^ per number of degrees 
of freedom (ndf ) are calculated based on the covariance matrix of the statistical errors. The 



results are compared with the OPAL data in Fig. 14. For details see text 
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Figure 1: Diagram of a photon-photon scattering process 
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Figure 2: The W^is distribution for all selected events at \^^^ = 183 GeV with W^is > 6 GeV 
after background subtraction. The data are compared to PHOJET (dashed line) and PYTHIA 
(dotted line). Only statistical errors are shown. The Monte Carlo background and e'j events 
which have been subtracted from the data are shown as the dashed-dotted histogram. 
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Figure 3: The selection efficiency defined by the ratio of the number of selected events, A'^sei, 
to the number of generated events, A'^gen, at a given generated invariant mass W for PHOJET 
(dashed line) and PYTHIA (dotted line) at v^g^ = 183 GeV. The lower curves give this ratio 
for the sum of the diffiactive events separately. 
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Figure 4: The relation between the visible hadronic invariant mass Wvis and the generated W for 
all selected PHOJET (dashed) and PYTHIA (dotted) Monte Carlo events at ^/s^^ = 183 GeV. 
The vertical bars show the standard deviation (spread) of the l^vis distribution in each bin. 
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Figure 5: The differential cross-section dcTcc/dW at (a) ^s^^ = 161 GeV, (b) \^^^ = 172 GeV 
and (c) \^^^ = 183 GeV for the anti-tagged two-photon events. Unfolding and acceptance 
corrections were done with PHO JET (dashed line) and with PYTHIA (dotted line). The 
statistical errors are smaller than the symbol size. 
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Figure 6: The A?7jnax distribution for all selected events at \/s^^ = 183 GeV with W^^is > 
6 GeV after subtracting the Monte Carlo background and e7 events. The A?7max distribution 
is compared to the full a) PHOJET and b) PYTHIA simulation. The diffractive component is 
shown separately. Only statistical errors are shown. 
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Figure 7: The distribution of the charged multiplicity rich for all selected events at i/i^g = 
183 GeV with VTvis > 6 GeV compared to PHOJET (dashed line) and PYTHIA (dotted line). 
The statistical errors are smaller than the symbol size. The subtracted Monte Carlo background 
and e7 events are shown as the dashed-dotted histogram. 
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Figure 8: The distribution of the thrust T for all selected events at \/Sgg = 183 GeV with 
H^vis > 6 GeV compared to PHOJET (dashed line) and PYTHIA (dotted line). The statistical 
errors are smaller than the symbol size. The subtracted Monte Carlo background and e'j events 
are shown as the dashed-dotted histogram. 
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Figure 9: The distribution of the energy Esw for all selected events at v^^g = 183 GeV with 
VTvis > 6 GeV and Egw > 1 GeV compared to PHOJET (dashed line) and PYTHIA (dotted 
line). The statistical errors are smaller than the symbol size. Due to the chosen scale the 
subtracted Monte Carlo background and e'j events are not visible. 
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Figure 10: The distribution of the energy Epi) for all selected events at v^ee ~ -'-^^ GeV with 
Wv,s > 6 GeV and Epo > 2 GeV compared to PHOJET (dashed line) and PYTHIA (dotted 
line). The statistical errors are smaller than the symbol size. Due to the chosen scale the 
subtracted Monte Carlo background and e7 events are not visible. 
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Figure 11: The distribution of the ratio PtZ-^vIs for all selected events at v^^,^ = 183 GeV 
with VTvis > 6 GeV compared to PHOJET (dashed line) and PYTHIA (dotted line). Only 
statistical errors are shown. The subtracted Monte Carlo background and e'j events are shown 
as the dashed-dotted histogram. 
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Figure 12: The distribution of the ratio Pl/E^is for all selected events at \^^^ = 161, 172 and 
183 GeV with W^is > 6 GeV compared to PHOJET (dashed line) and PYTHIA (dotted line). 
The statistical errors are smaller than the symbol size. Due to the chosen scale the subtracted 
Monte Carlo background and e'j events are not visible. 
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Figure 13: The total cross-section a^^(W) for the process 77 -^ hadrons. The OPAL 
measurement is compared to measurements by PLUTO 0, TPC/27 0, PEP/27 1, MDl 0] 
and L3 f^. The statistical and the systematic errors are added in quadrature. The normalisation 
uncertainty of 7% due to the extrapolation to Ql, Q2 = GeV^ is not included. The data are 
compared to a simple factorisation ansatz based on a Donnachie-LandshofF fit to total cross- 
sections [^] (solid line). The dashed-dotted line is the eikonalised mini-jet model by Godbole 
and Panchieri |^, the dotted line is the model of Schuler and Sjostrand [l^ and the dashed 
line is the model of Engel and Ranft M. 
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Figure 14: The total cross-section a^^{W) for the process 77 — * hadrons. Different Regge 
parametrisations have been fitted to the OPAL data. Only the statistical errors are shown. 
The vertical lines at the top of the figure delineate the bin boundaries. 
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